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Abstract 

Single particles and jets in deeply inelastic scattering at low x are measured with 
the HI detector in the region away from the current jet and towards the proton 
remnant, known as the forward region. Hadronic final state measurements in this 
region are expected to be particularly sensitive to QCD evolution effects. Jet cross- 
sections are presented as a function of Bjorken-x for forward jets produced with a 
polar angle to the proton direction, 0j et , in the range 7° < 6 iet < 20°. Azimuthal 
correlations are studied between the forward jet and the scattered lepton. Charged 
and neutral single particle production in the forward region are measured as a func- 
tion of Bjorken-x, in the range 5° < 8 < 25°, for particle transverse momenta larger 
than 1 GeV. QCD based Monte Carlo predictions and analytical calculations based 
on BFKL, CCFM and DGLAP evolution are compared to the data. Predictions 
based on the DGLAP approach fail to describe the data, except for those which 
allow for a resolved photon contribution. 
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1 Introduction 



The electron-proton collider HERA has made possible measurements of deeply inelastic 
scattering (DIS) in new kinematic regions: the regions of large four-momentum transfer 
Q 2 (up to Q 2 ~ 10 4 GeV 2 ) and small Bjorken-x (x ~ 10~ 4 ). A diagram illustrating 
the QCD evolution of a low x DIS event in the infinite momentum frame is shown in 
Fig. [TJ; a parton in the proton is able to undergo a QCD cascade resulting in several 
parton emissions before the final parton interacts with the virtual photon. These partons 
convert into hadrons which can then be detected and measured. The measurement of the 
hadronic final state between the proton remnant and the struck quark has been advocated 
over the past years as a technique for studying the dynamics of the QCD cascade in small 
x DIS 0,0,1. 

In this analysis data from the HI experiment are used to study single particles and 
jets in the forward region, defined as polar angle, 6, less than 25° in the laboratory 
reference frame measured with respect to the proton beam direction. A complication 
in using the hadronic final state to investigate QCD evolution arises due to the soft 
hadronisation processes. It has been suggested that this can be reduced by studying 
forward jet production |], ||, |!|. Such an analysis, based on a data sample containing 
approximately one tenth of the statistics used in this study, was reported in M . A similar 
analysis has been reported recently in . Single high-p r particle production as a test of 
the underlying QCD dynamics has also been proposed ||. A first analysis of transverse 
momentum (p T ) and pseudorapidity (rj) spectra was made using charged particles in 
the pseudorapidity region 0.5 < rj < 2.5 in the 7*7? centre of mass frame 0, accessing 
particles with angles down to 8° in the laboratory frame. The advantages of studying 
single particles, rather than jets lie in their independence of any choice of jet finding 
algorithm and also the potential to reach smaller angles than would be possible with jets 
with broad spatial extent. For single particles however, fragmentation effects are expected 
to be more significant. Thus both measurements are complementary. 

Several prescriptions for the QCD dynamics in the region of phase space towards the 
proton remnant have been proposed and are compared to the HI data in this study. These 
include QCD parton evolution schemes such as the classical DGLAP (Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi) || evolution, the low x specific BFKL (Balitsky-Fadin-Kuraev- 
Lipatov) || evolution equation and the CCFM (Ciafaloni-Catani-Fiorani-Marchesini) flOU 
evolution equation, which forms a bridge between the BFKL and DGLAP approaches. 
When the squared transverse momentum of a parton emitted in the ladder (Fig.0) is 
larger than Q 2 the resolved photon picture becomes applicable, in which case the photon 
can be considered to interact via its own resolved partonic component. 

2 QCD calculations and phenomenology 

DGLAP evolution, in which the relevant evolution parameter is InQ 2 , has been suc- 
cessfully tested in large Q 2 processes and provides an excellent description of the scale 
dependence of the structure function of the proton. The BFKL evolution equation, for 
which the relevant evolution parameter is hxl/x, describes scattering processes in QCD in 
the limit of large energies and fixed, but sufficiently large, momentum transfers. At lowest 



4 




Figure 1: Parton evolution in the ladder approximation. The kinematics of forward jets 
in DIS events are indicated. 

order the BFKL and DGLAP evolution equations effectively resum the leading logarith- 
mic contributions a s lnl/x and a s lnQ 2 respectively which, in an axial gauge, amounts 
to a resummation of ladder diagrams of the type shown in Fig. [Q. Since it is expected 
that at small enough x the lnl/x terms dominate the evolution, the BFKL equation 
should be applicable in this region. In the leading log DGLAP scheme the parton cascade 
follows a strong ordering in transverse momenta k\ n S> k\ n _ x 3> ... 3> feyi, while there 
is only a soft (kinematic) ordering for the fractional momenta x n < x n _\ < ... < x\. 
In the BFKL scheme, the cascade follows a strong ordering in fractional momentum 
x n <C x n _i <C ... <C X\, while there is no ordering in transverse momentum The 



transverse momenta perform a kind of random walk in kx space ; kn is close to kn-x, 
but it can be both larger or smaller ||12|| . For the CCFM evolution equation the relevant 



quantity is the angular ordering in the cascade. 

In a model that includes the resolved photon component in DIS the photon can be 
considered to interact via its own resolved partonic component. This internal photon 
structure leads to parton-parton scattering that can give rise to high E T jets, and emulate 
a random walk behaviour in kr ||13|| . This is usually theoretically treated by ascribing 



a structure function to the photon and allowing DGLAP evolution to take place from 
partons in both the photon and the proton. The ladder shown in Fig.|l] will then consist 
of two smaller ladders, one between the hard scattering and the photon, and one between 
the hard scattering and the proton. Consequently there is no longer a strict ordering in 
kx from the photon to the proton vertex, as obtained using the DGLAP picture without 
allowing virtual photon structure. 

Predictions for final state observables are available from Monte Carlo models, based 
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upon QCD phenomenology. LEPTO 6.5 [TJ|] incorporates first order QCD matrix ele- 
ments with matched leading-log DGLAP parton showers to simulate higher-order radia- 
tion (MEPS). For this model, the factorisation and renormalisation scales are set to Q 2 . 
The MEPS scheme fails to account for the di-jet rates when E^ et ^Q 2 Jl5], |16| in the kine- 



matic range studied in this paper, although this is not unexpected since for E^ et >Q 2 the 
parton from the proton is able to probe the structure of the photon. Using the concept of a 
resolved photon, as suggested above, the measured di-jet rates can be described fl3|, |i~6 



The Monte Carlo RAPGAP 2.06 |T7| can include such a resolved photon contribution 
in addition to the standard "direct" process. For the study presented here the SAS-2D 
parametrisation |T8[ of the virtual photon parton densities was used, and Q 2 + p?, with 



Pt the transverse momentum of the partons in the hard scattering, was used for the 
factorisation and renormalisation scales. 

ARIADNE 4.08 [0 provides an implementation of the Colour Dipole Model (CDM) 
of a chain of independently radiating dipoles formed by emitted gluons p0| . Unlike that 



in LEPTO, the parton cascade generated in the CDM is not kr ordered f2TJ[ . Since all 
radiation is assumed to come from the dipole formed by the struck quark and the remnant, 
photon-gluon fusion events have to be added and are taken from the QCD matrix elements. 



The linked dipole chain (LDC) model |22j is a reformulation of the CCFM equation and 
redefines the separation of initial and final state QCD emissions using the CDM. This 
overcomes technical problems in calculating non-Sudakov form factors for parton splitting 
functions. Calculations of the hadronic final state based on this approach are available 
within the LDCMC 1.0 ]23[ Monte Carlo which matches exact first order matrix elements 
with the LDC-prescribed initial and final state parton emissions. Both the LDCMC and 
ARIADNE use p\ as the factorisation and renormalisation scales. 

For the hadronisation of the partonic system, all of the above models use the Lund 
string model as implemented in JETSET p3]. All of these models are used with the 



GRV-HO [25] proton parton densities. QED corrections are determined with the Monte 



Carlo program DJANGO EH|, using the CDM for the final state. 



Alternatively, analytical BFKL calculations P7] , |28| for forward jet production are 
available at the parton level. These calculations are based on asymptotic expressions 
derived from the BFKL equation at leading order (LO). In order to compare to single 
particle spectra a different approach has been adopted. The normalization of the cal- 
culation was fixed by comparing the prediction for the forward jet cross sections with 



the data, as reported in [28|. Once this normalization was fixed, predictions for particle 



spectra were made for the kinematic range given below, using fragmentation functions. 

Also available are exact calculations of jet cross-sections at the parton level , at fixed 
order a 2 in the strong coupling constant, for direct or pointlike interactions of the pho- 
ton with the proton, as implemented in the DISENT [^] package. A comparison with 
DISENT for di-jet cross-sections has been reported in fT6|| . 

Background from photoproduction processes, those events for which Q 2 ~ in which 
the electron remains undetected in the beampipe and a fake electron is detected in the 
hadronic final state, has been studied using the PHO JET []30| Monte Carlo program. This 
generator contains LO QCD matrix elements for hard subprocesses, a parton showers 
model and a phenomenological description of soft processes. 
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3 The HI detector 



A detailed description of the HI apparatus can be found elsewhere |3l| . The following 



section briefly describes the components of the detector relevant to this analysis. 

The hadronic energy flow and the scattered electron are measured with a liquid 
argon (LAr) calorimeter and a backward electromagnetic lead-scintillator calorimeter 



(BEMC) respectively. The LAr calorimeter |32[] extends over the polar angle range 
4° < 9 < 154° with full azimuthal coverage. It consists of an electromagnetic section 
with lead absorbers and a hadronic section with steel absorbers. Both sections are highly 
segmented in the transverse and the longitudinal direction, in particular in the forward 
region of the detector, having about 45 000 channels in total. The total depth of both 
calorimeters varies between 4.5 and 8 interaction lengths in the region 4° < 9 < 128°. 
Test beam measurements of the LAr calorimeter modules show an energy resolution of 



<Je/E « 0.50/y-E [GeV] © 0.02 for charged pions The hadronic energy measure- 

ment is performed by applying a weighting technique in order to account for the non- 
compensating nature of the calorimeter. The absolute scale of the hadronic energy is 
presently known to 4%, as determined from studies of the transverse momentum (pr) 
balance in DIS events with forward jets [Q. 

The BEMC (depth 22 radiation lengths or 1 interaction length) covers the backward 
region of the detector, 151° < 9 < 176°. The primary task of the BEMC is to trigger on 
DIS processes with Q 2 values ranging from 5 to 100 GeV 2 , and to measure the scattered 
electron in DIS events. For scattered electrons with an energy larger than 11 GeV, as 
used in this analysis, the trigger efficiency is better than 99%. The BEMC energy scale 
for electrons is known to an accuracy of 1% |E5|. Its resolution is given by <Je/E = 
0.10/ \[E [GeV] © 0.39/£[GeV] © 0.017 gfj. 

The calorimeters are surrounded by a superconducting solenoid providing a uniform 
magnetic field of 1.15 T parallel to the beam axis in the tracking region. Charged particle 
tracks are measured in the central tracker (CT) covering the polar angular range 20° < 
9 < 160° and the forward tracking (FT) system, covering the polar angular range 5° < 
9 < 25°. The CT consists of inner and outer cylindrical jet chambers, z-drift chambers 
and proportional chambers. The jet chambers, mounted concentrically around the beam 
line, cover a range of polar angle of 20° < 9 < 165° and maximally provide 65 space 
points in the radial plane for tracks with sufficiently large transverse momentum. The 



achieved resolutions £1X6 CT pj, 

/Pt « 0.009-pt [GeV] ©0.015 and a e = 20 mrad 0. The FT 
provides measurements of tracks in the range 5° < 9 < 25°. This device has three identical 
sections, each containing a series of proportional, and radial and planar drift chambers, 
specifically arranged to facilitate triggering and track reconstruction at low polar angle. 
For tracks fitted to the interaction vertex the resolution in px has been demonstrated to 
be (t Pt /pt ~ 0.02 ■ p T [GeV] © 0.1 and the angular resolution, erg, to be better than 1 
mrad [P7j]. 

A backward proportional chamber (BPC), in front of the BEMC with an angular ac- 
ceptance of 155.5° < 9 < 174.5° serves to identify electron candidates and to precisely 
measure their direction. Using information from the BPC, the BEMC and the recon- 
structed event vertex the polar angle of the scattered electron is known to about 1 mrad. 

The luminosity is measured using the reaction ep — > ejry with two TICl/TlBr crystal 
calorimeters, installed in the HERA tunnel. The electron tagger is located at z = —33 m 
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and the photon tagger at z — —103 m from the interaction point in the direction of the 
outgoing electron beam. 



4 Data selection and corrections 



Experimental data for this analysis were collected by HI during the 1994 running period, 
in which HERA collided 27.5 GeV electrons^ with 820 GeV protons. Integrated lumi- 
nosities of the data samples used for the jets, n° and charged particle measurements are 
2.8, 2.0 and 1.2 pb _1 respectively. These differences arise due to the importance of select- 
ing running periods with optimal experimental conditions for each of the measurements 
presented here. 

DIS events are selected ]5o via the scattered electron which is experimentally defined 



as a high energy cluster, i.e. a localised energy deposit in the BEMC, with a cluster 
radius less than 5 cm and with an associated hit in the BPC The scattered electron is 
required to satisfy E e > 12 GeV and 9 min = 156° (for the particle analysis), or E e > 11 
GeV and 9 min = 160° (for the jet analysis), 9 min < 9 e < 173° where E e and 9 e are 
the energy and angle of the scattered electron. All results quoted here are corrected to 
their respective range. The kinematic variables are determined using information from the 
scattered electron: Q 2 = <iE e E' e cos 2 (6 l e /2) and y — 1 — (E e /E' e ) sin 2 (6 l e /2) where E e is the 
incident beam energy. The scaling variable Bjorken-x is related to these quantities via the 
square of the centre of mass energy s: x = Q 2 / {ys). Further reduction of photoproduction 
background and the removal of events in which a high energy photon is radiated off the 
incoming electron are achieved by requiring J2j (Ej — p z ,j) > 35 GeV with the sum 
extending over all detected particles j in the event, and by the requirement y > 0.1, 
respectively. The events are required to have an event vertex which is within 30 cm of 
the nominal vertex position. 

All distributions of data shown in this paper have been corrected bin-by-bin for de- 
tector effects, including geometrical acceptance, for QED radiative effects, and for the 
detection efficiencies of tt° mesons, charged tracks and jets respectively. The correc- 
tions were determined using events generated with the ARIADNE Monte Carlo program 
and a full simulation of the HI detector response. 



4.1 Forward 7r selection 

The 7r°-mesons are measured using the dominant decay channel tt° — > 27. The tt° candi- 
dates are selected in the region 5° < 6 m < 25°, where 6 n is the polar angle of the produced 
7r°. Candidates are required to have an energy of x 7T =E n /E p > 0.01, with E p the proton 
beam energy, and a transverse momentum with respect to the beam axis, pT n , greater than 
1 GeV. Here, pt v is taken to be equal to the transverse energy Et^ = E v sin^. At the 
high 7T° energies considered here, the two photons from the decay cannot be separated, but 
appear as one object (cluster) in the calorimetric response. Photon induced showers are 
selected by measuring the shower shapes of candidate clusters. The selection criteria are 
based on the compact nature of electromagnetic showers as opposed to those of hadronic 

1 In 1994 the incident lepton at HERA was a positron although we use the generic name 'electron' for 
the incident and scattered lepton throughout this paper. 



8 



origin. The very fine segmentation of the LAr calorimeter in the forward region (with a 
cell size of 3.5 x 3.5 cm 2 and four-fold longitudinal segmentation for the electromagnetic 
calorimeter) makes possible a detailed study of the transverse and longitudinal energy 
spread and the energy of each cluster. The high particle density in the forward direction 
leads to overlapping showers of electromagnetic and hadronic origin. Candidates of this 
type are largely rejected by the following requirements. 

A 7T° candidate is required to have more than 90% of its energy deposited in the 
electromagnetic part of the LAr calorimeter. A "hot" core consisting of the most energetic 
group of continuous electromagnetic calorimeter cells of a cluster, which must include the 
hottest cell, is defined for each candidate |3*9|. More than 50% of the cluster energy is 



required to be deposited in this core. The lateral spread of the shower is quantified in 
terms of lateral shower moments calculated relative to the shower's principal axis |39| and 



required to be less than 4 cm pO ]. The longitudinal shower shape is used as a selection 
criterion via the fraction of the shower's energy deposited in each layer of cells in the 
electromagnetic part of the calorimeter. The precise specifications of these layers can be 
found in |59|]. The value obtained by subtracting the fraction of energy deposited in the 
fourth layer of cells from that in the second layer as seen from the interaction point is 
required to be larger than 0.4, thereby selecting showers which start to develop close to the 
calorimeter surface and are well contained in the electromagnetic part of the calorimeter, 
as expected for showers of electromagnetic origin. With this selection 1673 7r° candidates 
are found in the kinematic range 5° < 6^ < 25°, x n > 0.01 and pxw > 1 GeV with a 
detection efficiency above 39%. 

Monte Carlo studies using a detailed simulation of the HI detector predict that less 
than 20% of the n° candidates are due to other misidentified hadrons. These calculations 
also show that less than 10% of the candidates originate from secondary scattering of 
charged hadrons with passive material in the forward region, where the amount of material 
between the interaction point and the calorimeter surface is largest. This contamination 
is subtracted bin-by-bin. Background from photoproduction has been studied with the 
PHOJET Monte Carlo model, following the procedure used in the measurement of F 2 



35], and is found to make no significant contribution to the distributions shown. 



4.2 Charged particle selection 

Charged particles are selected using the forward (FT) and central (CT) tracking chambers, 
in the region 5° < 9 cp < 25°, where 9 cp is the polar angle of the charged particle 
produced. Candidates are required to have an energy of x cp = E cp /E p > 0.01 and a 
transverse momentum with respect to the beam axis, pr cp , greater than 1 GeV. Basic 
quality criteria are used to select well measured tracks originating from the primary 
interaction point. Tracks measured in the CT are required to have at least 10 associated 
hits and a radial track length of over 10 cm. Selected tracks in the FT must have at 
least one planar segment confirmed by another segment reconstructed in either the planar 
or radial chambers. Furthermore, the impact parameter in the plane transverse to the 
beamline of the track from the primary vertex is required to be less than 5 cm when 
measured at the vertex position in z. Where a track has been found in both the FT 
and CT, a combined momentum measurement is used. For this measurement, as for the 
7r°-mesons, the high particle multiplicity in the forward direction constitutes the main 
challenge. It gives rise to a large flux of soft particles produced in secondary interactions 



9 



in passive material both in the FT and preceding it, leading to a degradation of detection 
efficiency and resolution. 

After the application of these cuts 2251 charged particles were selected for this study. 
The efficiency of reconstructing a central track and associating it to the primary event 
vertex is over 95% ETJ. The equivalent forward track efficiency is over 40% p7| for tracks 



well contained within the FT. For tracks produced at low values of 9 (9 < 10°), an 
efficiency of 25% is obtained. 

The proportion of selected charged particles that are produced in secondary interac- 
tions in the passive material in and around the HI detector is negligible for central tracks. 
For the forward tracks it is estimated by a detailed simulation of the HI detector to be 
maximally 10%. This contribution is subtracted bin-by-bin. The assumed material dis- 
tributions within the detector in the forward region have been verified by measuring the 
distribution of the rate of 7 — > e + e~ conversion pairs throughout the forward region of 
HI. 



4.3 Forward jet selection 

The region in which the forward jet measurements in this paper are made is chosen 
such that the phase space for jet production arising via DGLAP evolution is suppressed 
compared to that available in the case of BFKL evolution. This is achieved by requiring 
^Tjet ~ Q 2 1 where k T ^ et is the transverse momentum squared of the forward partonic jet. 
Due to the strong ordering in k T , this leaves little room for DGLAP evolution. In addition, 
the momentum fraction of the jet, x- ict = Ej et /E p , should be as large as possible, whereas 
the momentum fraction x of the quark struck by the virtual photon should be small 
in order to maximize the phase space for jet production from BFKL evolution which is 
governed by the ratio x- ]Ct /x. Here Ej et and E p are the energies of the jet and the incoming 
proton respectively. The rate of events at small x with a jet at high x- }et should be much 
higher for the BFKL than for the DGLAP scenario, as was shown in 



A simple cone algorithm [R3] is used to find jets in the HI LAr calorimeter, requiring an 



Et larger than 3.5 GeV in a cone of radius R = \J Arj' 2 + A<p' 2 = 1.0 in the space of pseudo- 
rapidity rj and azimuthal angle <fi in the HERA frame of reference. This jet algorithm was 
found to give the best correspondence between jets at the detector and hadron level for 
this type of analysis j42|. We have measured the cross-section for events which have a 



"forward" jet defined by x jet = E ict /E p > 0.035, 0.5 < p T]e JQ 2 < 2, 7° < 9 ]ct < 20° and 
Prjet > 3.5 and 5 GeV, where prjet is the transverse momentum of the jet, taken to be 
equal to Et- A total of 1945 events contained at least one jet satisfying these criteria. 

The transverse energy flow around the forward jet axis, averaged over all selected 
events, is shown versus At] and A0 with respect to the reconstructed jet axis in Fig. |2|, 
for various prjet ranges. Distinct jet profiles are observed which are reasonably described 
by the DJANGO model with full detector simulation. 

Backgrounds from photoproduction and radiative events (for which the current jet is 
boosted into the forward region) have been studied with Monte Carlo and with real data. 
The proportion of radiative events was reduced by requiring a cluster with energy larger 
than 0.5 GeV and pseudorapidity between —1.3 and 1.5 units of rapidity away from the 
forward jet axis, to tag the hadronic activity of the current quark. A fraction of the 
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Figure 2: The average transverse energy flow around the forward jet axis over all events 
is shown for two different minimum prjet values as a function of A<ft (left ), integrated over 
IA77I < 1 and as a function of A77, (right), integrated over |A0| < 1 rad, as shown by the 
insets in the top figures. Here At] and A0 are measured with respect to the reconstructed 
jet axis. Also shown is the DJANGO expectation with full HI detector simulation. 



background events can be experimentally measured (tagged) using the detectors of the 
luminosity system: for radiative events with a photon emitted collinearly with the incident 
electron, the photon can be detected with the photon tagger; for photoproduction events 
the electron which is scattered through a small angle can be detected in the electron 
tagger. Using data samples with tagged events it is possible to check the Monte Carlo 
based estimates of the photoproduction background. This is found to be less than 3% 
and subtracted from the data presented below. 



The forward jet analysis has also been repeated using an unfolding procedure f34jl 



based on the Bayes unfolding method |}43| . The results agree very well with the bin-to-bin 
corrected data. 



5 Results 

In Fig. single particle spectra are presented as a function of Bjorken-x for 7r°-mesons and 
charged particles. The 7r°-mesons are shown for 5° < 6 n < 25° in two bins of x„ (0.01 < x n 
< 0.015 and 0.015 < i T < 0.15) for pxw > 1 GeV and for the higher threshold of > 2 
GeV in the x n range 0.01 < x T < 0.15. All distributions are normalized to the number 
of DIS events, N, in the kinematic region specified in section ^ and with the additional 
requirement that Bjorken-x is in the range 0.0002 < x < 0.00235. The full errors are the 
quadratic sum of the statistical (inner error bars) and systematic uncertainties. For the 
7T° measurements the systematic errors contain the model dependence of the correction, 
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Figure 3: The single particle spectra in Bjorken-x are shown for charged particles and 
7r°-mesons produced in the polar angle range 5° < 6 < 25°. In the upper right plot 
ARIADNE with JETSET was used to calculate the contribution of charged pions to 
the charged particle measurement which makes possible, when divided by two, a direct 
consistency check of the two measurements. For comparison, four different Monte Carlo 
models are overlaid, as well as an analytical calculation labelled BFKL(KLM) based on 
p8f . n n is the number of 7r°-mesons and N is the number of DIS events that fall into the 
specified kinematic range. The full errors are the quadratic sum of the statistical (inner 
error bars) and systematic uncertainties. 

estimated from the difference in correction factors obtained using LEPTO and ARIADNE, 
which leads to errors of between 20-30%; a 4% variation of the electromagnetic energy 
scale in the forward region gives a 5% uncertainty; the variation of tt° selection cuts 
contributes 5-10% systematic error; and the variation of tt° acceptance also produces a 
5-10% uncertainty. 

The dominant source of systematic uncertainty for the charged particle measurements 
stems from the dependence of the FT efficiency on the forward particle multiplicity, which 
leads to a large systematic error due to the model dependence of the correction procedure. 
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This gives typical point-to-point errors of between 10% and 25% in the analysis presented 
here. These numbers were evaluated using simulated LEPTO and ARIADNE events. 
Discrepancies between the simulation and the FT detector response give rise to further 
point-to-point systematic errors of about 15%. Other sources of systematic effects, such 
as the BEMC energy scale uncertainty and photoproduction background are included in 
the errors but are small compared to those quoted above. 

The data in Fig. |3] clearly show that the production of 7r°-mesons with high px increases 
towards small Bjorken-x. The binning was chosen to match the x resolution. Only for 
the lower bin in i T and px is a measurement of charged particles available. To check 
the consistency of both measurements the charged particles results have been converted 
to charged pions spectra using the charged pion fraction predicted by ARIADNE with 
JETSET, which describes the charged particle data well. From isospin symmetry the 
neutral pion spectrum is expected to be half of the total charged pion spectrum. The 
resulting charged pion spectrum, divided by a factor two, is overlaid on Fig. [|b, and is 
found to be in good agreement with the ir° measurement. The model predictions shown 
will be discussed below. 

The forward jet data for pxjet > 3.5 GeV and pxjet > 5 GeV are shown as a function 
of Bjorken-x in the range 0.0001 < x < 0.004 in Fig. |j. In the case that there are two jets 
in the event which satisfy the jet selection criteria (approximately 1% of the forward jet 
events, see below), the jet with the largest px^ct was taken. It was shown that with the 
HERA kinematics and detector constraints, which determine the smallest reachable angle 
for #j et , the largest differences between BFKL and fixed order calculations are expected 



for low values of pxjet, namely around 3.5 GeV |44|]. Reducing this value further leads 
to excessively large hadronisation corrections as predicted by ARIADNE and dubious 
jet selections. However, to study the sensitivity of this study to increasing pxjet, cross- 
sections with prjet > 5 GeV have also been measured. The data points in Fig. f| present 
the cross-section as measured in the given x bin for the kinematic and jet cuts defined 
above. The bin purity, i.e. the fraction of events in a bin in x with a jet at the detector 
level which also have a jet at the hadron level in that bin, is between 25% and 50% for 
the lower Pxjet threshold value and somewhat larger for the higher Px)et threshold value. 

The data shown in Fig. § clearly rise with decreasing x, except for the smallest x bin, 
for which the rise is hidden due to the narrow kinematic region left with the present cuts. 
The total cross-section for forward jets with PTjet > 3.5 GeV in the defined kinematic 
range is 531 ± 17 (stat.)1ill(syst.) pb. 

The systematic errors on the forward jet data points shown in Fig. [| and in Table [1| 
are dominated by the 4% uncertainty of the hadronic energy scale, leading to shifts in the 
jet rates of 12%. Another important systematic contribution results from the dependence 
on the model used to correct the data and is up to 7% (the data shown are corrected with 
the ARIADNE model). The energy scale of the BEMC affects the cross-section by 3-4% 
in Figs. H and [5]. Other systematic effects have been found to be small compared to the 
above. 

The Monte Carlo models described above are compared to the data in Fig. |4j. The 
data rise faster with decreasing x than LDCMC and LEPTO, which are based on the 
CCFM equation and LO DGLAP dynamics respectively. The prediction of ARIADNE is 
slightly below the jet data and is slightly above the single particle data. ARIADNE does 
however reproduce the rise towards low x. Unlike LEPTO, this model does not include 
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Figure 4: The forward jet cross-section as a function of Bjorken-x for two Prjet cuts: 
3.5 GeV and 5.0 GeV. The errors shown are the statistical and systematic uncertainties 
added in quadrature. In (a) and (c) the curves are model calculations (full line RAPGAP, 
dashed line LDCMC, dotted line LEPTO and the dashed dotted line ARIADNE). In (b) 
(d) the full lines are analytic LO BFKL calculations at the parton level (no jet algorithm 
was applied), the dashed lines 0(ot 2 s ) calculations using DISENT with a cone algorithm 
applied. 

a kx ordered cascade and it has been suggested that these calculations will be close to 
the BFKL prediction pT]]. It has however also been hypothesised ||3 that the sizeable 



QCD radiation predicted by this model in the forward region stems not from the lack of a 
krp ordered cascade but from an unconventional enhancement of QCD-Compton processes 
at low x in which a hard gluon is radiated after the parton-photon interaction. The 
ARIADNE predictions also depend on the photon and target size extension parameters 
used in the calculations^. 

The RAPGAP model, which is based on LO DGLAP dynamics but has in addition a 



2 The default values PARA(10)=1.0 and PARA(14)=1.0 have been used. Changing these parameters 
to 1.5 (0.5) reduces (increases) the cross-sections by approximately 15-30% 
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Event Sample 


Forward jet cross-section [pb] 


HI 1994 Data 


531 ± 17 {stat.)t 8 i(sy st.) 


ARIADNE 4.08 


425 


RAPGAP 2.06 


491 


LEPTO 6.5 


259 


LDCMC 1.0 


262 



Table 1: Cross-sections for HI data and QCD calculations for inclusive forward jets 
satisfying x ]et > 0.035, 0.5 < p 2 T] JQ 2 < 2, 7° < jct < 20° and p Tjet > 3.5 GeV in the 
kinematic range E' e > 11 GeV, 160° < e < 173° and y > 0.1. 



contribution from resolved photons, gives a good description of both the single particle 
and jet data. The RAPGAP Monte Carlo allows the study of the scale dependence of 
the result. Sensible variations of the choice of scale usedf^in the leading order calculation 
were found to reduce the resolved photon cross-section calculations by up to 40%. The 
direct cross-section is more stable with respect to the choice of scale and changes by less 
than 10%. A comparison of the integrated forward jet cross-section for data and for these 
Monte Carlo models is given in Table |I[ The Monte Carlo cross-sections are based on 
event samples with numbers of events which are a factor 2 to 3 larger than those in the 
data samples. 

The forward jet data, shown in Fig. |] are also compared with recent analytical LO 
BFKL calculations an d a fixed order 0(a 2 ) jet calculation at the parton level. No 
jet algorithm is used for the BFKL calculation while a cone algorithm j|6) has been used 



for the fixed order calculation. The BFKL calculation is above the measurements for the 
data with both prjet thresholds of 3.5 and 5 GeV. The fixed order prediction for direct 
interactions of the photon with the proton is well below the data and very close to the so 
called "Born" calculation in JZ7[] , for which no gluon emissions are allowed in the BFKL 
ladder. The latter corresponds to the quark box diagram seen by the photon, and one 
single gluon at the low end of the diagram as is explained in f2"7fl . Variations of the fixed 
order cross section predictions of about 40% are found when the scale (Q 2 ) is varied by a 
factor of four. 

Hadronisation corrections and kinematic constraints (for the BFKL calculations) may 
change the predictions made at the parton level. The former are unfortunately strongly 
model dependent and at the lowest x values LEPTO and ARIADNE predict a change of 
the jet cross-section due to hadronisation by factors (number of jets at parton level/ num- 
ber of jets at hadron level) of 1.0 and 0.7 for PTjet > 3.5 GeV, and 1.5 and 1.1 for Prjct > 5 
GeV respectively. The kinematic constraints are expected to lower the predicted cross- 
sections |47j] . It should also be noted that the BFKL calculations contain a typical power 
growth of the cross-section with x with a power of about 0.5. Recently, next-to-leading 
order BFKL calculations for the inclusive cross-sections have become available |48| which 
demonstrate large corrections to this power, substantially reducing it. The question of 
how these large corrections will change the numerical BFKL predictions for the HERA 
forward jets cross-sections has not yet been addressed; in fact a recent paper [49| ar- 



gues that the next-to-leading order reduction of the BFKL exponent might be somewhat 



The scale was changed from + Q 2 to 4p 2 
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smaller than naively expected. Furthermore, the next-to-next-to-leading order corrections 
may also be large. In view of these uncertainties, the interpretation of possible deviations 
of experimental data from fixed order matrix element calculations as being due to BFKL 
effects has become more difficult. 

BFKL calculations folded with fragmentation functions are also available for the high- 
est momentum 7r° bins. This prediction is a factor of two above the data although the 
shape is in well described. This is shown in Fig. |3] in which the prediction, divided by 
two [pO[l , is compared to the data. Despite this normalisation discrepancy the calculations 
are seen to describe well the x dependence of the 7r° measurements. 




Figure 5: The azimuthal correlation A0 (in degrees) between the jet and the scattered 
electron, for two different x regions. The curves are model calculations (full line RAPGAP, 
dashed line LDCMC, dotted line LEPTO and the dashed dotted line ARIADNE). 

It has been suggested that a study of the azimuthal correlation A(f> between the scat- 
tered electron and the forward jet [[27]] may be a sensitive discriminant between the BFKL 
and DGLAP evolution schemes. A stronger decorrelation at low x is expected in the 
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Event Sample 


Forward di-jet cross-section (pb) 


HI 1994 Data 


6.0 ± 0.8(stat.)t^ii(syst.) 


ARIADNE 4.08 


6.7 


RAPGAP 2.06 


6.0 


LEPTO 6.5 


1.8 


LDCMC 1.0 


2.1 



Table 2: Cross-sections for HI data and different model calculations for forward di-jet 
events, each satisfying x jet > 0.035, 0.5 <p\- ]et /Q 2 < 2, 7° < # jct < 20° and p Tjct > 3.5 GeV 
in the kinematic range E e > 11 GeV , 160° <6 e < 173° and y > 0.1. 



BFKL case than in that of DGLAP due to the dilution of the original correlation result- 
ing from multi-gluon emission. The data are shown in Fig. [5] for the smaller and higher x 
regions. The data exhibit the expected levelling of the distribution when x decreases al- 
though all models follow this tendency, suggesting that this is not a very sensitive variable 
in this kinematic range. 

Of the 1945 forward jet events selected, 52 have a second jet which fulfill the forward 
jet criteria. This corresponds to a total cross-section of 6.0 ± 0.8(stat.)±3.2(syst.) pb, or 
roughly 1% of the forward jet cross-section. Recent LO BFKL numerical studies at the 
parton level [STJ predict that 3% of events will have two or more forward jets. We have 



also checked that this discrepancy between the data and the BFKL calculation remains if 
the ratio of Q 2 to Pxjet i s relaxed to 0.5 <p^ et /Q 2 < 5. However, the BFKL calculations, 
as before, include neither the effects of hadronisation nor those of jet algorithms. Ta- 
ble H compares the forward di-jet cross-section in data with the QCD models. ARIADNE 
and RAPGAP give di-jet rates in agreement with the data, while the other models pre- 
dict lower values and, owing to the current precision of the measurement, cannot yet be 
excluded. 



6 Summary 



Measurements of the cross sections for the production of forward jets and high trans- 
verse momentum single particles in low x DIS have been presented. Models whose dy- 
namics are based on different schemes for QCD evolution have been compared with the 
results. The conclusions using the forward jet and the single particle measurements are 
in agreement. 

Models implementing the traditional DGLAP evolution and including only direct pho- 
ton interactions grossly underestimate the amount of perturbative radiation required. Also 
a 0(a 2 ) jet calculation predicts less jet production. Models which include also the re- 
solved photon component successfully describe all the data presented. BFKL calculations 
at leading order describe the strong rise of the forward single jet and particle cross- sect ions 
at low x. While BFKL effects are expected to be prominent in this low x domain complete 
BFKL calculations have still to be performed at next-to-leading order before the data can 
be interpreted in this light. 
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Predictions of a model based on an implementation of CCFM evolution, which should 
smoothly interpolate between the DGLAP and BFKL regimes, give a poor description of 
all measured distributions. Calculations implementing the Colour Dipole Model provide 
sufficient QCD radiation to match the data. 
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